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Chaotic Hypothesis and Universal 


Large Deviations Properties 
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Abstract. Chaotic systems arise naturally in Statistical Mechanics 
and in Fluid Dynamics. A paradigm for their modelization are smooth 
hyperbolic systems. Are there consequences that can be drawn simply 
by assuming that a system is hyperbolic? here we present a few model 
independent general consequences which may have some relevance for the 
Physics of chaotic systems. 

Keywords and Phrases: Chaotic hypothesis, Anosov maps, Reversibility, 
Large deviations. Chaos 

§1. Chaotic motions.* 

A typical system exhibiting chaotic motions is a gas in a box whose particles 
interact via short range forces with a repulsive core, e.g. a hard core. No hope to 
ever be able to solve the evolution equations. 

In the very simple case of pure hard cores it has been possible to prove, 
mathematically at least in some cases, that the system is ergodic, [Si], [Sz], but 
ergodicity in itself is only a beginning of the qualitative theory of the motion. A 
similar situation arises in Fluid Mechanics: is a qualitative theory of Turbulence 
possible as, clearly, there are hopes to be able, in the near future, to prove an 
existence-uniqueness theorem but there is no hope for exact solutions of Navier 
Stokes equations? 

Equilibrium Statistical Mechanics is a brilliant example of a very successful 
quantitative theory derived from a comprehensive qualitative hypothesis, the er¬ 
godic hypothesis. The key to its success is a general expression for the probability 
distribution p, on phase space M providing us with the statistics of the motions 
corresponding to given values of the macroscopic parameters determining the state 
of the system. 

The statistics p is defined in terms of the time evolution map S via the 
relation: 



T-l 


( 1 . 1 ) 
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for all smooth observables F and for almost all, in the sense of volume measure 
on M, initial data x G M. In Equilibrium Statistical Mechanics the distribution 
fjL is identified with the uniform distribution on the surface of constant energy 
(the macroscopic state of the system being detemined by the volume V of the 
container box and by the energy U), which is an obviously invariant distribution 
by Liouville’s theorem of Hamiltonian Mechanics: this is a necessary consequence 
of the ergodic hypothesis. 

The success of Equilibrium Statistical Mechanics can be traced back to the fact 
that the ergodic hypothesis provides us with a concrete general, model independent, 
expression for the statistics of the motions. An expression that can be used to 
derive relations among time averages of various observables without even dreaming 
of ever being able to actually compute any of such averages. 

The Boltzmann’s heat theorem, the positivity of compressibility and specific 
heat are simple, but great, examples of such relations. They are relations which 
hold for any model, provided one makes the ergodicity hypothesis, see [Gal]. A 
classical argument that can be used to derive the heat theorem (i.e. the second law 
of Thermodynamics) from ergodicity is provided us by Boltzmann, see Appendix 
A2 and [Ga2]. 

Consider a mechanical system: viewing its phase space as a discrete set of 
points the ergodic hypothesis says that motion is a one cycle permutation of the 
points. Given a initial datum with energy U and with volume V we define temper¬ 
ature the time average of kinetic energy T = (AT) and pressure the time average of 
the derivative of the potential (/? with respect to the volume V (note that the force 
acting on the particles consists of the internal pair forces and of the force that the 
walls exercize upon the particles which depends on the position of the walls, hence 
it does change when the volume varies). Here and below (F) will denote the time 
average of the observable F. 

A general elementary property of a system whose motion on each energy 
surface is a single periodic motion is that if one calls p = {dv^p) then: 


dU -\-pdV 

-—-= exact 

T 


( 1 . 2 ) 


which means that if the energy U and a parameter V on which the potential 
depends (it will be the volume in our case) are varied by dU and dV respectively 
then the differential in (1-2) is exact. 

An elementary classical calculation shows that p, see Appendix A2, in the 
case of a gas in a box, has the meaning of average force exercized per unit surface 
on the walls of the container as a consequence of the particles collisions: thus we 
see that the ergodic hypothesis plus a general, trivial, identity among the averages 
of suitable mechanical quantities yields a relation (“equality of cross derivatives”) 
holding without free parameters. 

The reason why such relation is physically relevant for macroscopic systems 
is that the time necessary for the averages defining T,p to be reached within a 
good approximation by the finite time averages of K, dvP is not the unobservable 
recurrence time {i.e. the superastronomic time for the system to complete a single 
tour of the energy surface U) but it is a much shorter physically observable time 
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(whose theory is also due to Boltzmann being the essence of the Boltzmann’s 
equation) because the quantities K, dyf have an essentially constant value on 
the energy surface if the number of particles is large (so that the average of such 
observables “stabilizes” very rapidly compared to the recurrence times). 

To summarize: a simple hypothesis allows us to find the statistics of the 
motions of an equilibrium system: this implies simple parmeterless relations among 
averages of physically relevant quantities {i.e. = du^) which are observable 

in large systems because such quantities average very quickly compared to the 
recurrence times (being practically constant on the surface of given energy if the 
system is large). 

Thus a natural question arises: is there anything analogous in Non Equilib¬ 
rium Statistical Mechanics? and in developed Turbulence? 

The first problem is “what is the analogue of the uniform Liouville’s distri¬ 
bution?”. This is a really non trivial question that, once answered, will possibly 
allow us to try to find relations between time averages of mechanical quantities. 
The nontriviality is due to the fact that as soon as a system is out of equilibrium, 
i.e. nonconservative forces act upon it, dissipation is necessary in order to be able 
to reach a stationary state. But this means that any model used will be necessar¬ 
ily described by an evolution equation which will have a nonzero divergence: so 
that phase space will necessarily contract, in the average, and the statistics of the 
motion will be concentrated on a set of zero Lebesgue volume, see [Ru3]. 

Ruelle’s proposal in the early 1970’s was that one should regard such systems 
as hyperbolic so that there would be a unique stationary distribution describing 
the statistics of almost all initial data (chosen with the uniform distribution on 
phase space), [Rul]. 

This principle has been interpreted in [GC] as the following: 

Chaotic hypothesis: A chaotic mechanical system can be regarded for practical 
purposes as a topologically mixing Anosov system. 

This means that the closure of the attractor is a smooth surface on which 
the evolution is a Anosov system: of course assuming Axiom A instead of Anosov 
would be more natural, particularly in few degrees of freedom systems, [Rul]. 
However I prefer to formulate the hypothesis in terms of Anosov system as frac- 
tality of the closure of the attractor seems to be of little relevance in systems with 
large number of degrees of freedom occurring in Statistical Mechanics. 

The locution practical purposes is deliberately ambiguous as we know that 
even in Equilibrium Statistical Mechanics the corresponding ergodic hypothesis 
may fail while its consequences, at least some of them, will not (like the heat 
theorem in a free gas or in a harmonic chain). 

The above physical discussion serves as a quick motivation of the mathemat¬ 
ical question: are there general properties shared by mechanical systems that are 
transitive or mixing Anosov systems?. 

In the next sections I provide some affirmative answer in the class of time 
reversible Anosov maps and of weakly interacting chains of Anosov maps. Recall: 
a time reversal symmetry for a dynamical system (M, S) is any volume preserving 
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difFeomorphism / such that: 


1^ = 1, IS=S-^I (1.3) 

Examples in Hamiltonian mechanical systems are the velocity reversal, or the 
composition of the velocity reversal and the parity symmetry, or the composition 
of the velocity reversal, parity symmetry and charge conjugation symmetry. In 
general a time reversal may be a symmetry quite different from the naive one that 
can be imagined, see [BG]. 

Hamiltonian systems on which further anholonomic constraints are imposed 
via Gauss’ principle of least constraint often generate systems which show a time 
reversal symmetry, see Appendix Al, thus providing the simplest examples. 

§2. Time reversible dissipative Anosov systems. Fluctuation theorem. 

We now study a C°° , topologically mixing, Anosov system (M, S) on a compact 
manifold M. 

Let M be a d-dimensional, C°°, compact manifold and let S' be a C°°, mixing 
(transitive would sufhce) Anosov diffeomorphism, [AA], [Si]. If IT“, IT® denote the 
unstable or stable manifold at a; G M, we call IT"’'^, IT®’"^ the connected parts of 
IT“, IT® containing x and contained in the sphere with center x and radius S. Let 
du, ds be the dimensions of IT“, IT®: d = du + dg- We shall take S always smaller 
than the smallest curvature radius of IT“, IT® for x G M. Transitivity implies that 
IT“, IT® are dense in M for all a; G M. 

The map S can be regarded, locally near x, either as a map of M to M or of 
IT“ to IT_g,j,, or of IT® to Wg^. The Jacobian matrices of the ’’three” maps will be 
d X d, du X du and dg x dg matrices denoted respectively dS{x), 9S(a;)„, dS{x)s. 
The absolute values of the respective determinants will be denoted A(a:), A„(a;), 
As{x) and are Holder continuous functions, strictly positive (in fact A(a:) is C°°), 
[Si], [AA], [Ru4]. Likewise one can define the Jacobians of the n-th iterate of S; 
they are denoted by appending a label n to A, A„, A^ and are related to the latter 
by the differentiation chain rule: 


An{x) = A{S^x), Au,u{x) = AuiS^x), 

n (2-1) 

Ag^ni^^Ag ji{S^x^<^ A^(x) — Au^nix^ Ag^u{x')\ri{sA) 

3=0 

and Xn{x) = is the ratio of the sines of the angles a{S'^x) and a{x) 

between IT“ and IT® at the points S'”a; and x. Hence Xn{x) is bounded above 
and below in terms of a constant B > 0: B~^ < Xn{x) < 5, for all x (by the 
transversality of IT“ and IT®). 

We can define the forward and backward statistics or “SRB distributions” 
/i+, fj,- of the volume measure /io via the limits: 
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lim „ 
T^oo T 


1 ^ FiS^'^x) = f ^^±idy)Fiy) = m±(F) 


( 2 . 2 ) 


k=0 


which exist for all smooth functions F on M and for all but a set of zero volume 
of initial points x, see [Si]. 

Therefore it is the probability distribution /i.+ that is the statistics p, of the 
motions (almost surely with respect to the volume measure po on M), see (1-1): 
it plays the role of the Gibbs distribution, or microcanoncial ensemble, of equilib¬ 
rium Statistical Mechanics. Hence we are looking for general properties of p+, 
independent of the system considered, if possible. 

Let A(x) = I det9S'(a;)|; let p± be the forward and backward statistics of the 
volume measure pq {i.e. the SRB distributions for S and 


Definition: The system {M, S) is dissipative if: 


— p±{dx)logA^^{x) = rj^ > 0 

Jm 


(2.3) 


Remarks: 1) Existence of a time reversal symmetry I, see (1.3), implies ri_^_ = rj_ 
and A(a;) = A“^(/a;); furthermore /1E“ = kk/a, and the dimensions of the stable 
and unstable manifolds ds , du are equal: du = ds and d = d^ + dg is even. 

2) if Au{x), As{x) denote the absolute values of the Jacobian determinants of S as 
a map of kk“ to Wg^ and of Wf to Wg„,, then A„(a;) = As(/a;) 

3) If a system (M, S) is dissipative then the system (Af', S') with M' = M x M 
and S'{x,y) = {Sx, S~^y) provudes us with an example, setting I{x,y) = {y,x), 
of a dynamical system in the general class of “reversible” Anosov maps considered 
in §1. 

From now on only reversible dissipative Anosov dynamical systems {M, S) will 
be considered: it is for such systems that it will be possible to derive general model 
independent properties. 

Definition: The “dimensionless entropy production rate” or the “phase space con¬ 
traction rate” at x € M and over a time t is the function er{x): 


t/2-1 


(a;) = V logA ^(S'^x) = 3^1ogA.^ \x) 

71 . T f ^ 71 . T 


V+r 


i=-T/2 


d+'i 


(2.4) 


with Ar{x)'^= Hence (see (2.2)) it is, with po-probability 1; 


r/2-l 


1 T ) f 

(eT)+=.^lim -Ve^(S'^a;)= / p+{dy)sr{y) = 1 


(2.5) 


From the general theory of Anosov systems, [Si], it follows that the /im¬ 
probability that p = er{x) is in the interval [p — J, p + J] can be written as 
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maxqg[p _5 for some suitably chosen function C(p) and up to a factor 

bounded by with Q < B < +oo. This is a deep result of Sinai that holds 
because the statistics can be regarded as a Gibbs distribution and one can 
use the large deviation theory for such distributions: see Appendix A3 below for 
details. Then the following theorem holds, see [GC]: 

Fluctuation theorem: The “large deviation function” C,{p) is analytic in an interval 
i—p*,+p*) with p* > 1 and verifies the relation: 


C(p) -C(-p) 

PV+ 


\p\ < P* 


( 2 . 6 ) 


i.e. the odd part of C,{p) is in general linear and its slope is egual to the average 
entropy creation rate. 


What one really checks, see [Ga3], is the existence of p* > 1 such that the 
SRB distribution verifies: 


p — S < 


lim ^— log 

T^oo ry_|_T 


p +{{ er { x ) €[ p - S,p + S ]}) 
p +{{ er { x ) e -[p- (5,p + (5]}) 


<p + S 


for all p. Ip I < p*. 


(2.7) 


The above theorem was first informally proved in [GC] where its interest 
for nonequilibrium statistical mechanics was pointed out. The theorem can be 
regarded as a large deviation result for the probability distribution /x_|_. Although 
I think that the physical interest of the theorem far outweighs its mathematical 
aspects it is useful to see a formal proof. A proof is reproduced in Appendix A3 
below: it is taken from [Ga3]. 

The relation (2.6) has been tested numerically in several cases: it was in fact 
discovered in a numerical experiment, see [ECM2], and tested in other experiments, 
see [BGG], [BCL], [LLP]. Why does one need to test a theorem? the reason is that 
in concrete cases not only it is not known whether the system is Anosov but, in 
fact, it is usually clear that it is not, see [RT]. Hence the test is necessary to check 
the Chaotic Hypothesis which says that the failure of the Aoosov property should 
be irrelevant for “practical purposes”. 

Another interesting aspect, that cannot be treated here for limitations of 
time, of the above theorem is that it can be interpreted as an extension to non 
zero forcing (i.e. rj_^_ > 0) of the Green-Kubo relations: see [Ga6]. 


§3. Fluctuations in large systems. 


An important drawback of the above fluctuation theorem, besides the reversibility 
assumption which is not verified in many important cases, is that it can be prac¬ 
tically verified, for physical as well as mathematical reasons, only in (relatively) 
small systems. 

In fact the logarithm of the entropy creation rate distribution r C (p) is, usually, 
not only proportional to r, i.e. to the time interval over which the entropy creation 
fluctuation is observed, but also to the spatial extension of the system, i.e. to the 
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number of degrees of freedom; so that it is extremely unlikely that observing p in 
a large system one can see a value p which is appreciably different from 1 (note 
that the normalizing constant rj^ in (2.4) is so chosen that the average of p in the 
stationary state is 1). 

For this reason in macroscopic (or just “large”) systems the phase space con¬ 
traction rate is essentially constant (and its physical interpretation is of strength 
of the friction) much as the density is constant in gases at equilibrium. Therefore 
one can hope to see entropy creation rate fluctuations only if one can define a local 
entropy creation rate pvo (x) associated with a microscopic region Vq of space. 

I now discuss, heuristically, why one should expect that a local entropy creation 
rate can be defined, at least in some cases, and verifies a local version of the 
fluctuation law (2.6). This is discussed in a special example, see [Ga7], as in 
general one can doubt that a local version of the fluctuation law holds, see [BCL]. 

The special example that we select is the chain of weakly coupled Anosov 
maps, well studied in the literature, [PS]. The system has a translation invari¬ 
ant spatial structure, i.e. it is a chain (or a lattice) of weakly interacting chaotic 
(mixing Anosov) system. This can be described as follows. 

Let be a dynamical system whose phase space M' is a product of 

27V -I- 1 identical analytic manifolds Mq: M' = and S' : M' M' is a 

small perturbation of a product map x ... x Aq = Aq on M' . We assume that 
(7Vfo,S'o) is a mixing Anosov systems. The size TV (an integer) will be called the 
“spatial size” of the system. 

For x,y,z G Mq let F:.{x, z, y) be analytic and such that z ^ Fg{x, z, y) is a 
map, of Mq into itself, e-close to the identity and e-analytic for jej small enough. 
We suppose that, \i x = {x-n, ■ ■ ■, xn) G M': 


{S' x)i = Fe{xi-i,Xi,Xi+i) O SoXi (3.1) 

where a;-i-(jv-i-i) is identified with x^^n (*-e. we regard the chain as periodic); we call 
such a dynamical system a ehain of interacting Anosov maps coupled by nearest 
neighbors. It is a special example of the class of maps considered in [PS].^ 

It is difficult, maybe even impossible, to construct a (non trivial) reversible 
system of the above form: we therefore (see [Ga3]) consider the system {M,S) 

where M = M' x M' and define Sq'^= S^i x (Ao)“^ and S'^= S' x (5")“^, called 
hereafter the free evolution and the interacting evolution, respectively. So that the 
system can be considered as time reversible with a time reversal map I{ x, y) = 
{y, x). Note that the inverse map to (3.1) does not have the same form. The 
map S is, however, still in the class considered in [PS] because it can be written 
as S{x, y)i = {S{x, 11)11,S{x, y)^) with: 

S{x, y)ii =Fs{x^-i,Xi,Xi+i) o SoXi 
S{x, y)i2 =Ge,t{y) o S^^yi 


^ In the paper [PS] it is assumed that also Sq (hence Sq) is close to the identity, e.g. within 
e: such condition does not seem necessary for the purposes of the present paper, hence it 
will not be assumed. 
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where G has “short range”, z.e. |G£( y )i — G£( y')i| is of order if y and y'coin¬ 
cide on the sites j with \ j — i\ < k. By definition the system (M, S) is reversible”, 
i.e. the volume preserving diffeomorphism / verifies (1.3) above. 

Therefore the points of the phase space M will be {x,y) = 

{x-N,y-N ,. ■., xn, Un)' however, to simplify notations, we shall denote them by 
X = (x-AT ,... ,xn), with Xj denoting, of course, a pair of points in Mq. 

If e is small enough the interacting system will still be hyperbolic, i.e. for 
every point x it will be possible to define a stable and an unstable manifolds 
IT® , IT“ , [PS], so that the key notion of “Markov partition”, [Si], will make sense 
and it will allow us to transform, following the work [PS], the problem of studying 
the statistical properties of the dynamics of the system into an equivalent, but 
much more familiar, problem in equilibrium statistical mechanics of lattice spin 
systems interacting with short range forces. The reader will recognize below that 
this method is the natural extension to chains of the method used in Appendix 
A3 to study a single Anosov system. 

The main notion that we want to introduce for our chain is the notion of local 
entropy creation rate ??vb(a;), the entropy creation rate inside a fixed finite set 
Vo C [—A^, iV] of Anosov systems among the 2N -|- 1 composing the chain. 

Definition: Fixed a point x = (..., x^+i,...) consider the map (3.1) as 

dcf 

a map of Xy^ = {xj)j(zVo = {x-e,.. .,xe) into: 

2LVo *^(* ■ ■ ’ t—1 j Vo ’ ? ■ ■ ■) Vb (^■^) 

defined by (3.1) for i € We call “local entropy production rate” as¬ 

sociated with the “space like box” Vq = [—£, F\ at the phase space point x = 
{... ,Xi-i,Xi,Xi+i,...) the quantity riVoix) equal to minus the logarithm of the 
determinant of the 2{2£ -I- 1) x 2{2£ 1) Jacobian matrix of the map (3.3). 


Given a finite region Vq centered at the origin and a time interval 
To, let ? 7 _|_ denote the average density of entropy creation rate, i.e. 7]+ = 
limvo.To^oo VVoiS^x), then we set: 

. iTol 

P= —rr^ Y] VVoiS^x), V = VoxTo (3.4) 

ri+\V\ ^ 

' 3 = -iVo\ 

where yvoix) denotes the entropy creation rate in the region Vq- 

Calling TTvip) the probability distribution of p in the stationary state 
i.e. in the SRB distribution, and assuming that the system is a weakly coupled 
chain of Anosov systems I shall show, heristically, that: 


Proposition: It is irvip) = e^(^h^l+O(| 0 V|) \dV\ denotes the size of the 

boundary of the space-time region V and C,{p) is a function analyticinp S {—p*,p*) 
for some p* > 1. And: 


C(p) - Ci-p) 

PV+ 

Cip) =rC{p), 


1, IpI < p* 


V+ = rr]+ 


(3.5) 
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where r is the ratio between the total “volume” (2N + 1) of the system and the 
volume Vo, i.e. the “global” and “local” distributions are trivially related if appro¬ 
priately normalized. 

Note that this implies that if Vq is an interval of length L = |Vb| and if 
H = ITqI then the relative size of the error and of the leading term will be, for 
some length R, of order {L + H)R compared to order LH. Hence a relative error 
0{H~^ + L~^) is made by using simply ({p) to evaluate the logarithm of the 
probability of p as defined by (3.4)). 

The interest of the above statements lies in their independence on the total 
size 2N + 1 of the systems and the relevance of the above proposition for concrete 
applications should be clear. 

It means that the fluctuation theorem leads to observable consequences if one 
looks at the far more probable microscopic fluctuations of the local entropy creation 
rate. One can test the relation (3.5) in a small region Vq even when the system is 
very large: in such regions the entropy creation rate fluctuations will be frequent 
enough to be observable and carefully measurable. These fluctuations behave, 
therefore, just as ordinary density fluctuations at equilibrium: also the latter are 
not macroscopically observable but they are easily observable in small volumes. 

The key results for the analysis leading to the above proposition are the papers 
[GC], [Ga3] and, mainly, [PS]: the latter paper provides us with a deep analysis of 
chains of Anosov systems and it contains, I believe, all the ingredients necessary 
to make the analysis mathematically rigorous: however I do not attempt at a 
mathematical proof here. The analysis is presented in Appendix A4 below. 

Other types of fluctuation theorems (concerning non SRB distributions) had 
been previously found, see [ES]; extensions to stochastic systems have been recently 
discussed, see [Ku], [LSj. 

Acknowledgments: / have profited of stimulating discussions with F. Perroni, 
who also helped with numerical tests of the above ideas, with F. Bonetto and D. 
Ruelle. This work is part of the research program of the European Network on: 
“Stability and Universality in Classical Mechanics”, # ERBCHRXCT940460; par¬ 
tially supported also by Rutgers University and CNR-GNFM. 
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Appendix A1: The Gauss’ minimal constraint principle. 

Let vj(i, = 0, 2 = {ij, Xj} be a constraint and let R{i, x) be the con¬ 
straint reaction and x) the active force, see also Appendix A1 of [Ga3]. 

Consider all the possible accelerations a compatible with the constraints and 
a given initial state x,x. Then R is ideal or verifies the principle of minimal 
constraint if the actual accelerations + Ri) minimize the effort: 

A N 

— iPi-'nT'taif ^ ■6a^=0 (ATI) 

for all possible variations S a ^ compatible with the constraint p. Since all possible 
accelerations following i, x are such that ^) • = 0 we can 

write: 


£li — ixiitli — C(dx.<p{ij 2i) = Q 


with a such that -i^fix, x) = 0, i.e. : 


a = 




id± (py 


(A1.2) 


(A1.3) 


which is the analytic expression of the Gauss’ principle, see [LA]. 

Note that if the constraint is even in the i ^ then a is odd in the velocities: 
therefore if the constraint is imposed on a system with Hamiltonian H = K + V, 
with K quadratic in the velocities and V depending only on the positions, and 
if on the system act other purely positional forces (conservative or not) then the 
resulting equations of motion are reversible if time reversal is simply defined as 
velocity reversal. 

The gaussian principle has been somewhat overlooked in the Physics literature 
in Statistical Mechanics: its importance has been only recently brought again to 
the attention, see the review [HHP]. A notable, though ancient by now, exception 
is a paper of Gibbs, [Gi], which develops variational formulas which he relates to 
the Gauss principle of least constraint. 


Appendix A2. Heat theorem for monocyclic systems. Evaluation of 

THE AVERAGE {dyv) ■ 

Consider a 1-dimensional system with potential :p{x) such that |:/3^(a;)| > 0 for 
a;I > 0, > 0 and :p(x) -I- oo (in other words a 1-dimensional system 

in a confining potential). There is only one motion per energy value (up to a shift 
of the initial datum along its trajectory) and all motions are periodic so that the 
system is monocyclic. Assume also that the potential :p(x) depends on a parameter 

y. 

One defines state a motion with given energy E and given V. And: 

U = total energy of the system = K + ip 
T = time average of the kinetic energy K 
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V = the parameter on which (p is suposed to depend 
p = — time average of dyp 

A state is parameterized by U, V and if such parameters change by dU, dV respec¬ 
tively we define: 


dL =—pdV, dQ = dU + pdV (^2.1) 

then: 

Theorem (Helmholtz): the differential (dU + pdV)/T is exact. 

In fact let: 


r^+iU,V) _ rx+(u,v) _ 

S = 2 log / \/K{x; U, V)dx = 2 log / a/C/ — (p{x)d: 

Jx-{U,V) Jx-(U.V) 


x+iU,V) 

x-(U,V) 


(A2.2) 


( 5 S' is the logarithm of the action), so that: 


S = 


JidU-dvp{x)dV)^ 

J ^ Vk 


(A2.3) 


and, noting that ^^dt, we see that the time averages are given by inte¬ 

grating with respect to and dividing by the integral of We hnd therefore: 


dS = 


dU +pdV 

f 


(A2.4) 


Boltzmann saw that this was not a simple coincidence: his interesting (and 
healthy) view of the continuum (which he probably never really considered more 
than a convenient artifact, useful for computing quantities describing a discrete 
world) led him to think that in some sense monocyclicity was not a strong assump¬ 
tion. 

In general one can call monocyclic a system with the property that there is a 
curve £ x(£), parameterized by its curvilinear abscissa £, varying in an interval 

0 < £ < L{E), closed and such that x{£) covers all the positions compatible with 
the given energy E. 

Let X = x{£) be the parametric equations so that the energy conservation can 
be written: 


-I- ip{x{£)) = E (A2.5) 

then if we suppose that the potential energy p depends on a parameter V and if 
T is the average kinetic energy, p = —{dvp) it is, for some S: 

dS = p=-{dvp), T = (K) (A2.6) 
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where (•) denotes time average. 

The above can be applied to a gas in a box. Imagine the box containing the 
gas to be covered by a piston of section A and located to the right of the origin at 
distance L: so that V = AL. 

The microscopic model for the pistion will be a potential Tp{L — ^) if a; = 
(?) Vi C) &re the coordinates of a particle. The function ^(r) will vanish for r > tq, 
for some tq, and diverge to +00 at r = 0. Thus tq is the width of the layer near 
the piston where the force of the wall is felt by the particles that happen to roam 
there. 

Noting that the potential energy due to the walls is (/? = ~ 

that = A~^dL^ we must evaluate the time average of: 


dL^{x) =- ij) {A2.7) 

3 

As time evolves the particles with in the layer within tq of the wall will feel the 
force exercized by the wall and bounce back. Fixing the attention on one particle 
in the layer we see that it will contribute to the average of dL^(x) the amount: 


- 4—2 r 

total time Jf^ 


— ^j)dt 


(A2.8) 


if to is the first instant when the point j enters the layer and ti is the instante when 
the ^-compoent of the velocity vanishes “against the wall”. Since is 

the ^-component of the force, the integral is — 2 m\^j\ (by Newton’s law), provided 
> 0 of course. 

The number of such contributions to the average per unit time are therefore 
given by Pwaii A J_^^^2mv f{v) v dv if pmaii is the density (average) of the gas 
near the wall and f(v) is the fraction of particles with velocity between v and 
V + dv. Using the ergodic hypothesis {i.e. the microcanonical ensemble) and the 
equivalence of the ensembles to evaluate f{v) it follows that: 


p = {dv<p) = PwaiiP ^ (A2.9) 

where = ksT with T the absolute temperature and ks the Boltmann’s con¬ 
stant. That the (A2.9) yields the correct value of the pressure is well known, see 
[MP], in Classical Statistical Mechanics; in fact often it is even taken as microscopic 
definition of the pressure. 

Appendix A3. A proof of the fluctuation theorem. 


(A) Description of the SRB statistics. 

A set if is a rectangle with center x and axes A“, if: 

1) A", A® are two connected surface elements of kF“, Wf containing x. 

2) for any choice of ^ S A“ and 77 S A® the local manifolds and 1U“’‘^ intersect 

in one and only one point x4,r]) = n The point x{f,ri) will also be 

denoted ? x 77 . 
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3) the boundaries 9A“ and dA^ (regarding the latter sets as subsets of and 
W^) have zero surface area on and W^. 

4) E is the set of points A“ x A^. 

Note that any x' G E can be regarded as the center of E because there are 
A'“, A'^ both containing x' and such that A“ x A® = A'“ x A'®. Hence each E 
can be regarded as a rectangle centered at any x' G E (with suitable axes). See 
figure. 



The circle is a small neighborhood of ai; the first picture shows the axes; the intermediate picture 
shows the x operation and the third picture shows the rectangle E with the axes 

and the four marked points are the boundaries dA^ and dA^. The picture refers to a two 
dimensional case and the stable and unstable manifolds are drawn as fiat, i.e. the A’s are very 
small compared to the curvature of the manifolds. The center x is drawn in a central position, 
but it can be any other point of E provided A^ and A^ are correspondingly redefined. One 
should meditate on the symbolic nature of the drawing in the cases of higher dimension. 


The unstable boundary of a rectangle E will be the set duE = A" x i9A®; the 
stable boundary will be dgE = 9A“ x A®. The boundary dE is therefore dE = 
dgE U duE. The set of the interior points of E will be denoted A pavement 
of M will be a covering £ = {Ei,..., Ej^) of M by A/” rectangles with pairwise 
disjoint interiors. The stable (or unstable) boundary dgS (or du£) of £ is the union 
of the stable (or unstable) boundaries of the rectangles Ef du£ = yjjduEj and 

ds£ = UjdgEj. 

A pavement £ is called markovian if its stable boundary ds£ retracts on itself 
under the action of S and its unstable boundary retracts on itself under the action 
of S~^, [Si], [Bo], [Rul]; this means: 

Sds£ c ds£, S-^du£ C du£ (A3.1) 

Setting Mjji = 0, j,j' S {!,... , A/”}, if SE^ n Ej, = 0 and Mjji = 1 otherwise we 
call C the set of sequences j = (jfe)^_g^, jk G {!,... , A/”} such that = 1. 

The transitivity of the system (M, S) implies that M is transitive: i.e. there is a 
power of the matrix M with all entries positive. The space C will be called the 
space of the compatible symbolic sequences. If £ is a markovian pavement and d is 
small enough the map: 


OO 

X: j_gC ^ x= Pi S-'^Ej^ G M (A3.2) 


k——(yc) 


is continuous and 1 — 1 between the complement Mq C M of the set N = 
U'^_^S'^d£ and the complement Co C C of X~^(N). This map is called the 


Documenta Mathematica • Extra Volume ICM 1998 • 1-1000 



14 


Giovanni Gallavotti 


symbolic code of the points of M: it is a code that associates with each x ^ N a 
sequence of symbols j which are the labels of the rectangles of the pavement that 
are successively visited by the motion S^x. 

The symbolic code X transforms the action of S into the left shift d on C: 
SX{ j) = X(d j). A key result, [Si], is that it transforms the volume measure /io 
on M into a Gibbs distribution, [LR], [Ru2], JIq on C with formal Hamiltonian: 


— 1 OO 

H{j_)= ^ + + (A3.3) 

k— — (yc) k—0 

where, see (2.1): 


h-{j) =-'^ogAs{X{j)), h+{j) = log Au{X{j)), 

- - - (A3.4) 

hoii) = ~ logsina(X(j )) 

If F is Holder continuous on M the function F{j) = F{X{j )) can be repre¬ 
sented in terms of suitable functions 4>fc(j_fc,... ,jk) as: 


^ii) = |«'fe(j-fc,...,jfc)I < (,96 (A3.5) 

fe=i 

where (/? > 0, A > 0. In particular h± (and ho) enjoy the property (A3.5) {short 
range). 

If are the Gibbs states with formal Hamiltonians: 

OO OO 

^ h+{d^j_), Y. (^3.6) 

k——oo fe=— OO 

the distributions g,± on M, images of via the code X in (A3.2), will be the 
forward and backward statistics of the volume distribution yto (corresponding to 
JIq via the code X), [Si]. This means that: 

T—1 

hm i ^ F{S^^x) = [ ii±{dy)F{y) = ai±(F) (A3.7) 

for all smooth F and for /ip-almost all x S M. The distributions /i± are often called 
the SRB distributions, [ER]; the above statements and (A3.6),(A3.7) constitute the 
content of a well known theorem by Sinai, [Si]. 

An approximation theorem for /i+ can be given in terms of the coarse 
graining of M generated by the markovian pavement St = VIL-t S~^S.^ If 
= <A[^_rp and Xj_r^^,,,jj, is a point chosen in the coarse grain 

set so that its symbolic sequence is obtained by attaching to the right 


^ Where V denotes the operation which, given two pavements E,£' generates a new pavement 
£ s' ■. the rectangles of £ V simply consist of all the intersections iJ n of pairs of 
rectangles E e £ and e' e £^. 
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and to the left of j-r, ■ ■ ■ ,jT arbitrary compatible sequences depending only on 
the symbols j±T respectively. We define the distribution by setting: 


-:--i 


fiT.riF) = / HT,T{dx)F{x) = 

Jm 

t/2-1 

n A. (5'= a:) 




(H3.8) 


k=-T/2 


Then for all smooth F we have: lim 7 ’>T-/ 2 ,r^oo = ^J■+{F)■ Note that 

equation (A3.8) can also be written: 


^J'T,T{F) 










(A3.9) 


where j ° S C is the compatible sequence agreeing with j-r, ■ ■ ■ ,jT between —T 
and T {i.e. X{j^) = Xj_^^...jj, G continued outside as above. 


Notation: to simplify the notations we shall write, when T is regarded as having a 
fixed value, q for the elements q = {j-T, ■ ■ ■ jjr) of {1,..., and Eg will 

denote and Xg the above point of Eg . 


Remark: Note that the weights in (A3.9) depend on the special choices of the 
centers Xg {i.e. of j°); but if a;^ varies in Eg the weight of Xq changes by at 
most a factor, bounded above by some B < oo and below by B~^, for all T > 0, 
and essentially depending only on the symbols corresponding to the sites close to 
±T. 


The last formula shows that the forward statistics of /io can be regarded as a 
Gibbs state for a short range one dimensional spin chain with a hard core interac¬ 
tion. The spin at k is the value of jfc G {!,..., A/”}; the short range refers to the fact 
that the function h+{j) = logA„(X(j)), (A„(x) being Holder continuous), can 
be represented as in (A3.5) where the play the role of’’many spins” interaction 
potentials and the hard core refers to the fact that the only spin configurations j 
allowed are those with = 1 for all integers k. 

(B) A Legendre transform. 

First the function (2.4) is converted to a function on the spin configurations 
J GC: 


eAj_)=er{X{j_)) = ^ ^ L{Aj_) 

k=-T/2 


(A3.10) 
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where L{j) = ^ log A^^(X( j)) has a short range representation of the type 
(A3.5). 

The SRB distribution is regarded (see above) as a Gibbs state with 
short range potential on the space C of the compatible symbolic sequences, associ¬ 
ated with a Markov partition £, [Si], [Ru2]. Therefore, by general large deviations 
properties of short range Ising systems ([La], [El], [01], there is a function ({s) real 
analytic in s for s G {—p*,P*) for a suitable p* > 0, strictly convex and such that 
lip <p* and [p — S,p + S] C {—p*,p*) we have: 

-log-p+{{er{j) £[p-S,p + S]})^.^^ max C(s) (A3.11) 

T ~ se |p—<5,p-l-<5] 

and the difference between the r.h.s. and the l.h.s. tends to 0 bounded by Dt~^ 
for a suitable constant D. The function C(s) is the Legendre transform of the 
function A(/3) defined as: 


m 


lim —log J ^_i_(d j ) 


(A3.12) 


i.e. A(/3) = max 5 g(_p. p.)(/3s -f C('S)), where the quantity p* can be taken p* = 
lim/ 3 ^+oo gP~^X{P) and the function A(/3) is a real analytic, [CO], strictly convex 
function of /3 G (—oo, oo) and /3~^A(/3) i.e. it is asymptotically linear. 

The above (A3.11) is a ’’large deviations theorem” for one dimensional spin 
chains with short range interactions, [La]. 

Hence it will be sufficient to prove the following; if Ip^s = [p — P + iJ]: 


1 M+({£r( j ) € Ip,s^T]{T)}) f < p + 6 + 77'(r) 

V+T °®71+({er(j) G/_p,5±p(^)}) \ > p - S -p'It) 


(A3.13) 


with p{t), rj'{ t) 


(C) Thermodynamic formalism informations. 


In this section X will denote a lattice interval, i.e. a set of consecutive integers 
X = (cc, a; -I- 1,..., a: -I- n — 1): hence it should not be confused with the code X of 
(A3.2). 

Let = {jx,jx+i,---,jx+n-i) if-A = (x,a;-b l,...,a:-bn-1) and n is 
odd, and call X = x + {n — l)/2 the center of X. If j G G is an infinite spin 
configuration we also denote j ^ the set of the spins with labels x € X. The left 
shift of the interval X will be denoted by i?; i.e. by the same symbol of the left 
shift of a (infinite) spin configuration j . 

Let = l^'^Hjx,jx+i,---,jx+n-i), and = 

h^+''{jx,jx+i,---,jx+n-i) be translation invariant, z.e. functions such that 
l^x{j) = lx{j) and h~^{j), and such that the functions h+{j), 

see (2.4), and L{j), see (A3.10), can be written for suitably chosen constants 
bi,b2,b,b': 
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= ^+(z) = X! 

17=0 17=0 (H3.14) 

\lx{j_^)\ \h+{j_^)\<be-^'- 

Then T£r{3_) can be written as Y.x^[-t/ 2 ,r/ 2 - 1 ] ^^(Ix^' 

Hence T£r{j ) can be approximated by Te^{j ) = Yl^^'^lxij where 
means summation over the sets X C [—^t+M], while 17 is in [—^r, ^r—1]. 
The approximation is described by: 

\T£f{i) - T£r{l)\ < 636 "'’^^ (H3.15) 

for suitable^ 63,64 and for all M >0. Therefore if Ip^s = [p — + S\ and M = 0 

we have: 


P+{{^t{x) G ^p,s}) 


— M+({^t € /p^^+bg/r}) 

— P+({^r € Ip^S—bs/r}) 


(A3.16) 


It follows from the general theory of 1-dimensional Gibbs distributions, [Ru2], 
that the 7 i+“probability of a spin configuration which coincides with j ^ ^ 

in the interval [—^r, \t\,^ is: 


=-E* 




E. 


-[-x/2,t/2] L J 


^^i[-r/2,T/2]^ 


(H3.17) 


where denotes summation over all the X C [— t/2,t/2 — 1]; the denominator 
is just the sum of terms like the numerator, evaluated at a generic (compatible) 
spin configuration j |_^^2 r/ 2 ] ’ ^ verifies the bound, [Ru2]: 


Bl^<Pill-r/2.r/2l^<B^ (^ 3 . 18 ) 

with Bi a suitable constant independent of j j _^/2 t/ 2 ] ^ explic¬ 

itly estimated in terms of b,b'). Therefore from (A3.16) and (A3.17) we deduce 
for any T > r/ 2 : 


€ Ip^s}) ^ /I+({E ^ ^p,S+b3/T}) ^ 

< B2PT,t{{^ G Ip^S+bs/A) < B2PT,t{{£t G -fp,5+2&3/r}) 


(A3.19) 


for some constant H 2 > 0 ; and likewise a lower bound is obtained by replacing B 2 
by BA and 63 by - 63 . 


One can check from (A3.14), that the constants 63,64 can be expressed as simple functions 


of fci, 62 . 

i.e. the spin configurations j' such that = jx, x G [—^r]. 
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Then ii p < p* and Ip^s C {—p*,p*) the set of the rectangles E S V-t 
with center x such that er{x) € Ip^s is not empty, as it follows from the strict 
convexity and the asymptotic linearity of the function A(/3) in (A3.12). 

We immediately deduce the lemma: 

Lemma 1: the distributions /i+ and T >^t, verify: 


1 A^+({£t(2 :) S /p,5=p2b3/T}) 

T^+ ^ IJ,+ {{£t{x) e -/p,5±2&3/r}) 


< 


log 


> - 


logS^ 




-L- log 
rri+ “ Ut 

log 7 


ut,t (l^T e Ip, 5 1) 


1 

rv+ 


{{St£ — Ip,s}) 
Wr,T (l^Te1) 


TP 


-£—Ip,s}) 

(A3.20) 


for Ip s C [—p*,p*] and for r so large that p + 5 + 2 & 3 /T < p*. 

Hence (A3.13) will follow if we can prove: 

Lemma 2: there is a constant b such that the approximate SRB distribution fj,T,T 
verifies: 


1 /i-T,T({gr € Ip,s}) ( < p + S + b/r 
rj+T ^ fJ,T,T{{£T &-Ip,5}) \>p-S-b/T 

for T large enough (so that 6 + b/T < p* — p) and for all T > t/2. 


(A3.21) 


The latter lemma will be proved in §4 and it is the only statement that does 
not follow from the already existing literature. 

(D) Time reversal symmetry implications 

The relation (A3.20) holds for any choice of the Markov partition £. Note 
that if f is a Markov pavement so is i£ (because iS = S~^i and iWtf = kkj^); 
furthermore if £i and £2 are Markov pavements then £ = £\\l £2 is also markovian. 
Therefore: 


Lemma 3: there exists a time reversal Markov pavement £, i.e. a Markov pavement 
such that £ = i£. 


This can be seen by taking any Markov pavement £0 and setting £ = £q\I i£o. 
Alternatively one could construct the Markov pavement in such a way that it 
verifies automatically the symmetry [G2]. Since the center of a rectangle Eg G £t 
can be taken to be any point Xq in the rectangle Eg we can and shall suppose 
that the centers of the rectangles in £t have been so chosen that the center of iEg 
is iXq , i.e. the time reversal of the center Xg of Eg . 

For T large enough the set of configurations q = j such that Sr (x) G 

Ip^s for all X G Eg is not empty® and the ratio in (A3.21) can be written, if is 
the center of Eg G £t, as: 


® Note that p* = sup^, limsup^_^_l_gQ sriS'''^‘^x) and let p e (—p* + i5,p* — 5); furthermore 
({s) is smooth, hence > — 00 , for all |s| < p*. 
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^e^{xg)elp,s ^u,Ti^q) ^e-rixg)elp,s q) 

'^e-rix q)G-Ip,S ^U,Ti^q) (s , ) 6/p ,5 ^ g) 


(H3.22) 


Define As^p(x) as in (A3.8) with s replacing u: then the time reversal sym¬ 
metry implies that Au,t{x) = Ag p{ix), see remark 2) following definition (B), §2.^ 
This permits us to change (A3.22) into: 


Xe.^ (x 
Xet (a: 

where the maxima are evaluated as q varies with erixq) G Ip^s. 

By (2.1) we can replace A^^^(x)A^ l(x) with A^ ^{x)B^^, see (A3.8), (2.4); thus 
noting that by definition of the set of q ’s in the maximum in (A3.23) we have 

log Ap ^{xq ) G we see that (A3.21) follows with b = ^ \ogB. 

Corollary: the above analysis gives us a concrete bound on the speed at which the 
limits in (2.6) are approached. Namely the error has order 

This is so because the limit (A3.11) is reached at speed furthermore the 

regularity of A(s) in (A3.11) and the size of rj{T),ri'{T) and the error term in 
(A3.21) have all order 

The above analysis proves a large deviation result for the probability distribu¬ 
tion /i+: since pL+ is a Gibbs distribution, see (A3.6), various other large deviations 
theorems hold for it, [DV], [El], [01], but unlike the above they are not related to 
the time reversal symmetry. 

Appendix A4: Heuristic proof of the local fluctuation theorem. 


--1 / 


,)eip,s ( < maxq A^l{x q)Ag ^{x q) 


,)&ip 


^s.rixq) [ > min, A„ p(a;q)A^ ^(a;,) 


(A3.23) 


(A) Markov partitions and symbolic dynamics for the chain. 

The reduction of the dynamical nonequilibrium problem of a weakly interact¬ 
ing chain of Anosov maps, see §3, to a short range lattice spin system equilibrium 
problem is the content of (A), (B) of this appendix, see [Ga7]. This is an extension 
of the corresponding analysis in Appendix A3 for the case of a single Anosov map: 
it is necessary to discuss it again in order to exploit the short range nature of the 
coupling and its weakness in order to obtain results independent on the size N of 
the chain. 

Let Vo = {Ei,...,E^^) be a Markov partition, see [Si], for the unper¬ 
turbed “single site” system (Mq x Mo,So x Sq ^). Then Vo^~^^ = {Ea}, 


Here it is essential that Au,t{x) is the expansion of the unstable manifold between the initial 
point S~'''Cx and the final point S’’’Cx-. i.e. it is a trajectory of time length t, which at 
its central time is in a;. 
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a = (p-AT,... ,pAr) with Ea = X X ... X E°^ is a Markov parti¬ 

tion of 

The perturbation, if small enough, will deform the partition into a 

Markov partition V for (M, S) changing only “slightly” the partition Vq . The 
work [PS] shows that the above “e small enough” mean that e has to be chosen 
small but that it can be chosen N -independent, as we shall always suppose in what 
follows. 

Under such circumstances we can establish a correspondence between points of 
M that have the same “symbolic history” (or “symbolic dynamics”) along 
under Sq and along V under S'; we shall denote it by h; see [PS]. 

_2A^+1 

The Markov partition Vq for So associates with each point x = 
{x-N, ■ ■ ■ ,xn) a sequence (cTij), i G [—N,N],j G (—00,00) of symbols so that 
i<^i,j)j^-oo is the free symbolic dynamics of the point Xi. We call the first label 
i of aij a “space-label” and the second a “time-label”. Not all sequences can 
arise as histories of points: however (by the definition of h, see above) precisely 
the same sequences arise as histories of points along Vq under the free evolution 
Sq or along V under the interacting evolution S. 

The map h is Holder continuous and “short ranged”: 

\h{x)i-h{x')^\<CY,e\^-^\-^'\xj-x)y (H4.1) 

3 

for some 7,7', G > 0 , [PS], if \x — y\ denotes the distance in Mq x Mq {i.e. in the 
single site phase space). 

Furthermore the code x <—> a associating with x its “history” or “symbolic 
dynamics” ct (x) along the partition V under the map S is such that, fixed j: 

2 L{-£)i = \i — j\ < £ ^ \xj — x'j\< Ce'^^ (H4.2) 

The inverse code associating with a history g_ a point with such history will be 
denoted x{g_). 

If X = (x-AT,... ,xn) is coded into ct(x) = ..., a^) = (crij), with 

i = —N,.. .,N, and j G (—00, -l-oo), the short range property holds also in the 
time direction. This means that, fixed iq: 

aij = a'i j ior \i - io\ < k,\j\ < p ^ \ x{q_)ip - x{a')ig\ < 

(H4.3) 

for some k,^,C > 0 , see lemma 1 of [PS]. The constants k,^,C,C',B,B' > 0 
above and below should not be thought to be the same even when denoted by the 
same symbol: however they could be a posteriori fixed so that to equal symbols 
correspond equal values. 

By construction the codes x <—> commute with time evolution. 

The sequences (ctij ) which arise as symbolic dynamics along Vo under the free 
single site evolution of a point Xi are subject to constraints, that we call “vertical”, 
imposing that r°. = 1 for all j, if denotes the “compatibility matrix” 
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of the “free single site evolution” (z.e. = 1 if the S'o x S'g ^ image of E^r 

intersects the interior of Ea' and =0 otherwise). We call the latter condition 
a “compatibility condition” for the spins in the t-th column. 

The mixing property of the free evolution immediately implies that a large 
enough power of the compatibility matrix has all entries positive. This means 
that for each symbol cr we can find semiinfinite sequences: 

(Jb{(j) a-i,ao = a), T° =1, for all z < 0 

n (AAA) 

crT((j) =(cr = (70,0-1 ,...), = 1, for all z > 0 

and defines two functions aB,crT, called “compatible extensions”, defined on the 
set {1,..., A/q} of labels of the single site Markov partition Vo, with values in the 
compatible semiinfinite sequences. 

In fact there are (uncountably) many ways of performing such compatible 
extensions “from the bottom” and “from the top” of the symbol a into semiinfinite 
compatible sequences of symbols. We imagine to select one pair aB, cft arbitrarily, 
once and for all, and call such a selection a “choice of boundary conditions” or “of 
extensions”, on symbolic dynamics, for reasons that should become clear shortly. 
All this seems unavoidable and it is closely parallel to the corresponding discussion 
in the analysis of the simpler case of a single Anosov system discussed in Appendix 
A3, see the discussion preceding (A3. 8 ). 

We shall therefore be able to “extend in a standard way” any finite compatible 
block® Q of spins: 

g_Q = {crij)i(zL,j(iK, L = {a - £,a +tj, K = {h - m,h +m) (A4.5) 

by setting = aB{cri,b-n)b-n-j iov j < b - n and atj = aT{(Ti,b+n)j-b-n for 
j > b + n. Here a, 6 , i, m are integers. 

In the free evolution there are no “horizontal” compatibility constraints; hence 
it is always possible to extend the finite block CTq = {aij)i^L,jGK to a “full 
spin configuration” sequence ((Tij)jg[_jv,Ar],je(-00,00)) obtained by continuing the 
columns in the just described standard way, using the boundary extensions cts, ar, 
above the top and below the bottom, into a biinfinite sequence and also by extend¬ 
ing the spin configuration to the right and to the left to a sequence with spatial 
labels running over the full spatial range [—N, A^]. One simply defines aij for i ^ L 
as any (but prefixed once and for all) compatible biinfinite sequence of symbols 
(the same for each column). 

The allowed symbolic dynamics sequences for the free dynamics (on Vo) and 
for the interacting dynamics (on V) coincide because the free and the interacting 
dynamics are conjugated by the map h, [PS]. Therefore the above operations make 
sense also if the sequences are regarded as symbolic sequences of the interacting 
dynamics, as we shall do from now on. 

To conclude: given a “block” CTq of symbols, with space-time labels (z,j) G 
Q = L X K, we can associate with it a point x G M whose symbolic dynamics 

® A block j), (i,j) e Q, is naturally said to be “compatible” if ^ o’* j+i “ ^ 

(b j) e Q such that (i,j + 1) is also in Q. 
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is the above described standard extension of a q. The latter depends only on the 
values of aij for j at the top or at the bottom of Q and, of course, on the boundary 
conditions as, <Jt chosen to begin with. 

(B) Expansion and contraction rates. 

Consider the rates of variation of the phase space volume, Ao(x), or, respec¬ 
tively, of the surface elements of the stable and unstable manifolds As(x) and 
Xu{x) at the point x: they are the logarithms of the Jacobian determinants 
dS{x), i9(a)S'(x), a = s,u, where 9 (q,) denotes the Jacobian of S' as a map of 
to Wg^ where a = u,s distinguishes the unstable manifold of x or the stable 
manifold of x: 

Aa(ir) = — log I det 9 (q)S( X )|, q; = 0,u, s (A4.6) 

where d(^o)S{x)'^= dS{x). 

A hard technical problem is to represent Xa{x) in terms of the “symbolic 
history” of x along V, i.e. in terms of compatible sequences a = (cTij) with 
i € {—N,N), j G (—cx), oo). The rates Aq-(x) can be expressed as: 

Aa(ai) =-log|det|^|j^„(^^ = (^4.7) 

— Lcl-N,N] 

where L is an interval in [—iV, TV] (with ±(A^ + 1) identified with =fA^)) [PS]. 

For a = 0 this can be done by noting that the matrix J = ^ has an al¬ 
most diagonal structure: J{x) = Jo{x){l + A(x)) where Jo{x) is the Jaco¬ 
bian matrix of the free motion Jq{x) = Jo(x-n) x Jo^X-n+i) x ... x Joix^) if 
X = {x-n,...,xn) and if D = ( ^et Jo(a:j)): 

det J = D • log(l+A( 2 =)) ^ ^ 

which leads to (A4.7) if one uses that the matrix elements Ap^g = 
Jo^{x)dxj,dx^J{x) are essentially local, i.e. bounded by i? for some 

7 ,G,B > 0 (see (3.1),(3.2), (A4.3)). 

For a = u, s (A4.7) can be derived in a similar way using also that: 

(1) the stable and unstable manifolds of x consist of points y which have 
eventually, respectively towards the future or towards the past, the same history 
of X, 

(2) they are described in a local system of coordinates around x = 
{... ,X-i,xo,xi,...) by smooth “short range” functions. Suppose, in fact, 
that on each factor Mq we introduce a local system of coordinates (a, /3) around 
the point Xi G Mq, such that the unperturbed stable and unstable manifolds are 
described locally by graphs {a,fs{a)) or (/„(/?), /J). 

The unperturbed stable and unstable manifolds will be smooth graphs 
(cci, fs(a.i)) or {fuiPi),Pi) with ai varying close to cti and Pi close to Pi, with 
(ai,Pi) being the coordinates of Xi. 
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Fixed a point x = {x-n, ■ ■ ■ ,xn) with coordinates {cti, I3^)i--N,...,N the per¬ 
turbed manifolds of the point x will be described by smooth (at least and in 
fact of any prefixed smoothness if e is sufficiently small) functions W^{a), fF“( (3 ) 
of a = {ai)i=-N,N or of P = {Pi)i--N,N which are “local”; i.e. if a and a' agree 
on the sites i — £,i + £ or if P and P' agree on the sites i — £p + £ then: 

\\W-{P)i-fu{Pi)\\c-<Ce, \\W-{PP-W^{P%\\c^<Ce^ 

- - - ^ (H4.9) 

\\W^{a)^ - fsiaP\\c2 < Ce, \\W^{aP - fF*(a')*llc= < Ce^ 

for some C > 0, see [PS] lemmata 1,2. Here the norms in the first column are 
the norms in as functions of the arguments P or respectively a , while the 
norms in the second column are norms evaluated (of course) after identifying 
the arguments of P (or a) and P' (or a') with labels j such that \i — j\ < £. 

(3) If we consider the dependence of the planes tangent to the stable and unstable 
manifolds W'^. , 1F“ at x we find that they are Holder continuous as functions of 
x : 

\{dWl)i-{dWl)i\<cY,e\^-^\^\x,-yP\ a = u,s (H4.10) 

3 

where {dW^)i denoted the components relative to the i-th coordinate of x of the 
tangent plane to and C, k, 7 > 0 . 

The above properties and the Holder continuity (A4.1), (A4.2), (A4.3) imply 
that the “horizontal potentials” 5^j^\xiP in (A4.7) are “short ranged”: 

a = u,s (A4.11) 

for some H, C, 7 > 0; we denote \L\ the number of points in the set L. 

We shall use the symbolic representation of x € M to express the rates 
A^“^(x). For this purpose let x = {xpi^-N^N and suppose that such x cor¬ 
responds to the symbolic dynamics sequence a = {g_j)fp_oo where = 
■ ■ ■ ,crN,j)- We denote the sequence = (cTj.jOieij=-oo,oo- 
Then aj^ does not determine Xj^ (unless there is no interaction, i.e. e = 0): 
however the short range property, (A4.3), of the symbolic codes and of the map 
h conjugating the free evolution and the interacting evolution shows that, if V 
is a larger interval containing L and centered around L, then the sequence g_j^, 
determines each point of Xj^ within an approximation < l-l^'l) 7 ^ Hence 

we can define ^“^(ct^,) so that: 

S‘t\xL) = 

L'ZPL 

L 

for some This leads to expressing Xa{x) in terms of the symbolic dy¬ 

namics of X and of the “space-localized” potentials ?^“^(ct^). 


l4“^(aJ|<H'(CV)l^l-i 


(A4.12) 
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Let Qn = L X K where K = [—n, n] is a “time-interval” and set 

(^ 4 . 13 ) 

if n > 1 and [ctq„] denotes a standard extension (in the sense of §3) of or 

just set Cq^ 'y\[n-Qo]) n = 0. We dehne >Cq(ctq) for Q = ^ x ^ and 
K not centered {i.e. K = {a — n, a + n), a 7 ^ 0) so that it is translation invariant 
with respect to space time translations (of course the horizontal translation invari¬ 
ance is already implied by the above definitions and the corresponding translation 
invariance of 

The remarkable property, consequence of the Holder continuity of the functions 
in (A4.6) and of the (A4.3),(A4.12), see [PS], is that for some j,k,B,C > 0: 

\C2i{aQ)\<BiCe^ye-^^ (A4.14) 

if i,j are the horizontal and vertical dimensions of Q. 

In this way we dehne a “space-time local potential” which is, by 

construction, translation invariant and such that, if A denotes the box A = 
[—N, A^] X [—M, M] the following representations for the rates in (A4.6) hold: 

-log|deta(„)5'2“+^(5'"^x)| = ^/:§(CTQ)-bO(|5A|) (A4.15) 

QcA 

where 0(|9A|) is a “boundary correction” due to the fact that in (A4.15) one 
should really extend the sum over the Q’s centered at height < M and contained 
in the inhnite strip [—N, A^j x [— 00 , 00 ] rather than restricting Q to the region A. 
Hence the remainder in (A4.15) can, in principle, be explicitly written, in terms 
of the potentials in the boundary term form usual in Statistical Mechanics 

of the 2-dimensional short range Ising model and it can be estimated to be of 
0(|i9A|) by using (A4.14). 

(C) Symmetries. SRB states and fluctuations. 

Besides the obvious translation invariance symmetry the dynamical system 
has a time reversal symmetry, this is the diffeomorphism I, see (1.3), which anti¬ 
commutes with S and Sy. 


IS = S-^I, ISo = Sor\ 1^ = 1 (A4.16) 

We can suppose that the Markov partition is time reversible, i.e. to each element 
of the partition V one can associate an element Ea-' = IE,y which is also an 
element of the partition. Here we simply use the invariance of the Markov partition 
property under maps that either commute or anticommute with the evolution S'. 
hence it is not restrictive, see [Ga5],[Ga3], to suppose that for each g_ one can 
dehne a ct' so that E^' = IE ^ . We shall denote such a' as I a or also — a. For 
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e = 0, z.e. for vanishing perturbation, the map / will act independently on each 
column of spins of a . This property remains valid for small perturbations; hence: 

= {a'ij} = {-CTi.j} (H4.17) 

i.e. time reversal simply reflects the spin configuration corresponding to a phase 
space point and changes “sign” of each spin. 

The functions Aq(x ) and their “potentials” -Cq(o;q) verify, as a consequence, 
if Q = X [—fc, k] is a centered rectangle: 

K{Ix) =-K’{x), ^gigiq) =-Cgi-a^g) (H4.18) 

where a' = s if a = u and a' = u if a = s, oi = 0 if a = 0. The above symmetries 
will be translated into remarkable properties of the SRB distribution. 

The “local entropy production rate” associated with the “space like box” Vq = 
[—(,(] at the phase space point x = {... ,Xi-i,Xi,Xi+i,.. .) has been defined in 
§3 in therms of the Jacobian matrix of the map S. We can likewise consider the 
corresponding Jacobian determinants of the restriction of the map S to the stable 
and unstable manifolds of x. Such determinants will depend not only from Xi, 
i € Vq, and on the nearest neighbors variables x±£ but also on the other ones Xk 
with |fc| > t' + 1: however their dependence from the variables with labels \k\ > £ 
is exponentially damped as by (A4.14). Thus we can define riy^,iiiy^ in a 

way completely analogous to py^ in (3.3). 

If we look at the average phase space variation rates rfy ^, lyy ^, between the 
time —£) and i? we can find, via a power expansion like the one in (A4.8) along the 
lines leading from (A4.8) to (A4.15), a mathematical expression as: 

VvoiiS-) — ''^*^q{3Lq) (A4.19) 

Q 

where the Q runs over rectangles Q centered at 0-time Q = [a—£, a+£] x [—fc, k] 
with [a — £,a -\- £] C Vq. This could be taken as an alternative definition of r]y^, as 
it is a rather natural expression. For our purposes, if H = Vq x [—'d,'!?], one needs 
to note that (A4.19) holds at least in the sense that: 


1 

Ro • {2d + 1) 




j=-i} 


1 

Ro • (2i? + 1) 


Qcv 


Oi\dV\) 

1^1 


(A4.20) 


i.e. expression (A4.19) can be used to compute the average local entropy creation 
rate in the space-time region V up to boundary corrections 0{\dV\) (that can be 
neglected for the purposes of the following discussion). 

We now study the SRB distribution /i: denoting by {F)+ the average value 
with respect to yi of the observable F we can say, see [Si], [PS], that if A = 
[-N,N] X [-T,T]: 


{F)+ 


lim 

T —»-oo 




(A4.21) 


Documenta Mathematica • Extra Volume ICM 1998 • 1-1000 



26 


Giovanni Gallavotti 


We want to study the properties of the fluctuations of: 

'+ Qcv Qcv 

for which we expect a distribution of the form 7ry(p) = const The 

SRB distribution gives to the event that p is in the interval dp the probability 
Trv{p)dp with: 


T^v{p) = const ^ (Ad.23) 

at fixed p 

and (defining implicitly C/“): 

E^q(^q)= E^q(^q)+ E J^U^Q)+Oi\dV\n-^)‘‘^^ 

QCA Qcv QCA/V (A4.24) 

— UY{ay) + + 0(\dV\ K 

with K > 0, having used the “short range” properties (A4.14) of the potential. 

In the sums in (A4.21) we would like to sum over ay and over g_j^iy as if 
such spins were independent labels. This is not possible because of the vertical 
compatibility constraints. However the mixing property supposed on the free 
evolution implies that the compatibility matrix T° raised to a large power R has 
positive entries. Hence if we leave a gap of width R above and below V we can 
regard as independent labels the labels aij with i in the space part Vq of the region 
V = Vo X [—I?,!?] and with \j\ > d + i?, by a distance > R above or below the 

de "f 

region V. Denoted V + R = Vq x [—d — i?, d + i?] remark that: 

E = UyiS-v) + UA/iv+R)i3LA/(v+R)) + 0{\dV\ {R + K ^)) (A4.25) 

QcA 


Hence, proceeding as in [GCl], we change the sum over (the dummy label) g_ in 
the denominator to a sum over — a^ and using C^,{— g_Q) = —Eq(ctq): 


7rv(p) 

TTv{-p) 


E 


at fixed p 


J^QCV — p^A/(V + R)^2.A/(V + R)') 


c,J2qcv “-^q(2:q) I/a/(v+b)((-2;^)a/(v+R)) 


E 


gO(|3v|) (^4 26 ) 


at fixed p 


with the summation being over the spin configurations in the “whole space-time” 
A, subject to the specified constraint of having the same value for p, i.e. the 
same average local entropy creation rate in the space-time region V. The latter 
expression becomes, since the labels a, — a^ (respectively in the numerator and 
denominator of (A4.26)) are independent dummy labels: 


Ea. + R)) 
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so that by the (A4.20), (A4.22) and since the symmetry relations above im¬ 
ply the relation Eqcy(^q(^q) +^q(5 : q)) = V rj+p, up to corrections of size 
0{\dV\K~^) we find, (note the repetition of the comparison argument given in 
[GC]): 

(A4.28) 

Ty(-p) 

yielding a local fluctuation law, i.e. the first of (3.5). The second line of (3.5) is a 
(simple) consequence of the above analysis but we do not discuss it here. 
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